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Abstract 



The Poyakov-Nambu-Jona-Lasinio (PNJL) model was developed recently, which includes both the chiral dynamics and 
(de)confinement effect and gives a good description of lattice QCD data. In this study we use the PNJL model to describe the 
- - - quark phase, and first use it to study the evolution of proto-neutron star (PNS) with a hadron-quark phase transition. Along the 
line of a PNS evolution, we take several snapshots of PNS profiles, presenting the fractions of different species, the equations of 
state (EOS), and the mass-radius relations at different stages. The calculation shows the mixed phase may exist during the whole 
evolving process, and the onset density of quark phase decreases with the radiation of neutrinos in the heating stage. In the cooling 

Onstage, the EOS of the mixed phase softens and the center density increases. In this process a part of nuclear matter transforms to 

^ quark matter, which may lead to a PNS collapsing into a black hole. 
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Since the matter in the core of compact stars are compressed 
to densities of several times of saturated nuclear density, it is 
expected that new degrees of freedom will appear in the inte- 
rior of these objects The 



used to describe cold neutron star without 113011 and with a color 
superconductivity structure 13 If . In 13211 . Fischer et al. dis- 



hadron-quark phase transition is one of the most concerned top- 
ics in modern physics related to heavy-ion collision experiment 
and compact star Because of the complication of full calcula- 
tion of Quantum Chromodynamics (QCD) and the lack of suf- 
ficient knowledge about the nonpertubtive and (de)confinement 
effect, it is difficult to apply full QCD calculation to describe 
the phase transition in astrophysics. Therefore, in literatures 
the hadron-quark phase transition related to compact star are 
usually described with the simplified, phenomenological MIT 
bag model or the effective chiral Nambu-Jona-Lasinio (NJL) 

taodei (e.g., iSiBIiilElEIBESEl). 

The NJL model with chiral dynamics is a prominent one in 
the application of astrophysics, but it lacks the confinement 
mechanism, one essential characteristic of QCD. Recently, an 
improved version of the NJL model coupled to Polyakov loop 
(PNJL) has been proposed ll20ll . The PNJL model includes both 
the chiral dynamic and (de)confinement effec t, g ivin g a goo d 



interpretation of lattice QCD data lH miHlll, E?!] 



We have recently studied the hadron-quark phase transition rel- 
evant to heavy-ion collision in the Hadron-PNJL model 1128 1. 
The calculation shows the color (de)confinement effect is very 
important for the hadron-quark phase transition at finite density 
and temperature, and improves greatly the results derived from 
the Hadron-N JL model 12911 . The PNJL model with a chemical 
potential dependent Poyakov effective potential has also been 
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cussed the evolution of core collapse supernova with the PNJL 
model (with diquark interaction and isoscalar vector interac- 
tion) and the possibility of the onset of deconfinement in core 
collapse supernova simulation. In their study selected proton- 
to-baryon ratio Yp and the Maxwell construction are taken for 
the PNJL model, which gives a relatively narrow mixed phase 
than that with Gibbs construction. In 13311 . the authors studied 
the possibility to probe the QCD critical endpoint during the 
dynamical black hole formation from a gravitational collapse. 
Several two-flavor quark models with different parameters are 
used in their study, and the calculation shows the Critical Point 
location of QCD has a strong dependence on quark models and 
parameters. A generalization with s quark is needed to get more 
reliable results for further investigation. 

In this Letter we will firstly use the newly improved PNJL 
model to describe the evolution of proto-neutron star from 
its birth with trapped neutrinos to neutrino-free cold neutron 
star (NS). The Gibbs criteria will be used to determine the 
mixed phase under the isotropic constraint with and without 
trapped neutrinos. The emphases are put on the evolution of 
proto-neutron star with a hadron-quark phase transition, the 
particle distributions along baryon density, the EOSs and mass- 
radius relations, as well as the star stability in different snap- 
shots during the evolution. 

A PNS forms after the gravitational collapse of the core of 
massive star with the explosion of a supernova. At the be- 
ginning of the birth of a PNS, the entropy per baryon is about 
one (5^1) and the number of leptons per baryon with trapped 
neutrino is approximate 0.4 (Yie = i'e + i'v,, - 0.4). In the fol- 
lowing 10-20 seconds, neutrinos escape from the star. With 
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the decrease of electron neutrino population, the star matter is 
heated by the diffusing neutrinos, and the corresponding en- 
tropy density increases, reaching to 5 - 2 when Yy^ ^ 0. Fol- 
lowing the heating, the star begins cooling by radiating neutrino 



pairs of all flavors, and finally a cold neutrons forms 1134 1351] . 

Along the line of a PNS evolution to the formation of a cold 
neutron star, we take several snapshots to study how the star 
evolves, especially with the appearance of quark degrees of 
freedom. The snapshots are taken with the following conditions 
(5 = 1, Fi, = + Fv. = 0.4), (5 = 1.5, i-L, = + Fv, = 0.3), 
(S = 2, K = 0) and (S = 0, Fv, = 0), similar with that 
used in 1341 13511 . At each snapshot of a PNS evolution, we 



take an isentropic approximation with which the temperature 
has a radial gradient in the star. There are also studies related to 
the prop erties of PNS based on isothermal approximation (e.g., 

EoHB )• 

For the star matter, the hadronic and quark phase are de- 
scribed by the non-linear Walecka model and the PNJL model, 
respectively. In the mixed phase between pure hadronic and 
quark matter, the two phases are connected to each other by 
the Gibbs conditions deduced from thermal, chemical and me- 
chanical equilibriums. For the hadron phase we use the La- 
grangian given in ll36ll in which the interactions between nu- 
cleons are mediated by cr, co, p mesons, and the parameter set 
GMl is used in the calculation. The details can be found in 
Refs. Sll. 

For the quark phase, we take recently developed three-flavor 
PNJL model with the Lagrangian density 
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where q denotes the quark fields with three flavors, m, d, and 
s, and three colors; m^) - diag(m„, m^, nis) in flavor space; 
G and K are the four-point and six-point interacting constants, 
respectively. The four-point interaction term in the Lagrangian 
keeps the S t/v(3)x5 t/A(3)xt/v(l)xt/A(l) symmetry, while the 
't Hooft six-point interaction term breaks the Ua{1) symmetry. 

The covariant derivative in the Lagrangian is defined as 
D/t - - iAfj. The gluon background field = ^'^Aq is 
supposed to be homogeneous and static, with Aq = g^J, Y' 
where ^ is SU{3) color generators. The effective potential 
'ZY(<1)[A], 3)[A], T) is expressed in terms of the traced Polyakov 
loop d) = (Tt,L)/Nc and its conjugate O = (Tr,.L' )/A^c- The 
Polyakov loop L is a matrix in color space 



L(x) = f exp 
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where /3 - l/T is the inverse of temperature and A4 - iAq. 
Different effective potentials were adopted in literatures 121 

3' 



37, 38 



^^ 



[]. The modified chemical dependent one 



asT^ ln(l - 6(D- + 80^ - 30)'*) 



(3) 



was used in 1130113 111 which is a simplification of 



(aoT^ + ciif/ + a2T-fr)^(t) 
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because the difference between O and <1> is smaller at finite 
chemical potential and <1) = O at yu = 0. In the calculation we 
will use the later one. The related parameters, oq - -1 .85, aj - 
-1.44 X lO-\ fl2 = -0.08, fl3 = -0.4, are still taken from [30|l , 
which can reproduce well the data obtained in lattice QCD cal- 
culation. 

In the mean field approximation, quarks can be taken as free 
quasiparticles with constituent masses M,, and the dynamical 
quark masses (gap equations) are obtained as 



Mi = Mi - 4G(f>i + 2K(/)j(f)k (i + i + k). 



(5) 



where tpi stands for quark condensate. The thermodynamic po- 
tential of the PNJL model in the mean field level can be derived 



Q = 
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where A(<D, O, - T) = 1 + 3<t>e-*^'-^''>/^ + 3(I)e-2(£-/J,)/r + 
g-3(£,-;:,)/r and A(€), O, + yL/,-, T) = 1 + 3€)e-<^'+''''/^ + 

The values of 0„, (pd, <ps, and O are determined by minimiz- 
ing the thermodynamical potential 
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All the thermodynamic quantities relevant to the bulk prop- 
erties of quark matter can be obtained from D.. Particularly, 
the pressure and entropy density can be derived with P - 
-(Q(r,;u) - Q(0, 0)) and S = -dQ/dT, respectively. 

As an effective model, the (P)NJL model is not renormaliz- 
able, so a cut-off A is implemented in 3-momentum space for 
divergent integrations. The model parameters: A - 603.2 Me V, 
GA2 = L835, KA^ = 12.36, ot„,^ = 5.5 and m, = 140.7 MeV, 
determined by fitting f„, M^, rriK and to their experimental 
values 1391], are used in the calculation. 

The Gibbs criteria is usually implemented for a complicated 
system with more than one conservation charge. The Gibbs 
conditions for the mixed phase of hadron-quark phase transition 
in compact star are 



jQ 



(8) 
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where ju,, are usually chosen with //„ and yu^. Under the (3 equi- 
librium with trapped neutrino, the chemical potential of other 
particles including all baryons, quarks, and leptons can be de- 
rived by 

l^i ^ bijUn - qil^e + qil^v,, (9) 

where bj and qi are the baryon number and electric charge 
number of particle species respectively. For the matter with 
trapped electron neutrinos, Yj^ = (7^ + Fy,,) - 0, we do not 
need to consider the contribution from muon and muon neu- 



trino 1134 13511 . For the neutrino-free matter (//y , = 0), both elec- 
trons and muons are included in the calculation. 

The baryon number density and energy density in the mixed 
phase are composed of two parts with the following combina- 
tions 

(10) 



P = (1 -X)Pb +Xpt 



and 



(11) 



where is the volume fraction of quark matter And the electric 
neutrality is fulfilled globally with 

qtotai = (1 -x) 2 q'pi +x^q'pi = 0. (12) 

i=B,l i=q,l 
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the smaller YpjYn will excite more neutrons to occupy higher 
energy levels, leading to a stiffer equation of state. 

Another point we stress is that the onset density of quark 
phase decreases with the escape of trapped neutrinos, i.e., 
trapped neutrinos delay the hadron-quark phase transition to a 
higher density. For the neutrino-trapped matter, the fraction of 
neutrino is enhanced with the appearance of quarks, which also 
affects the neutrino opacity. For neutrino-free cases, the lep- 
ton population in cold neutron star matter (5 = 0, Yy^ - 0) is 
smaller than that of PNS {S - 2, Fy,, = 0), especially after the 
appearance of quarks at high density. 

The largest center densities of proto-neutron stars at the first 
three snapshots taken above are - 5.23, 5.00, and 4.54 po^ 
respectively. The center density decreases with the radiation 
of neutrinos, because the star expands when the inner matter is 
heated by the escaped neutrinos. When neutrinos are free, the 
PNS begins cooling by radiating neutrino pairs of all flavor and 
the star shrinks until the formation of cold neutron star. During 
the cooling stage, the center density of the star increases. The 
center density of cold neutron star is p^, = 4.88 p^. The varia- 
tion of the center densities during a PNS evolution is easier to 
understand by combining the mass-radius relations that will be 
given latterly. 
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Figure 2: EOSs of PNS and NS matter at several stages of the star evolution. 
The dots mark the range of the mixed phase. 



Figure 1 : Relative fractions of dilferent species as functions of baryon density at 
several snapshots of a PNS evolution. The upper (lower) panels are the results 
with (without) trapped neutrinos 

In Fig.[T] we display the relative fractions of different species 
as functions of baryon density in several snapshots along the 
evolution of a proto-neutron star, from (5=1, Y^e - Y^ + Yy^ = 
0.4), (S = 1.5, Y^e ^Ye + Yy^ = 0.3), (5 = 2, Yy^ = 0) to cold 
neutron star with (S = 0, Yy^ - 0). Comparing the upper pan- 
els with trapped neutrino with the lower panels without trapped 
neutrino, we find that the fraction of trapped neutrino affects the 
proton-neutron ratio Yp/Y„ at lower density before the appear- 
ance of quarks. For the hot PNS matter, Yp/Y„ with rich is 
larger than that with poor . According to the Pauli principle. 



We present the equations of state of the star matter at differ- 
ent evolving stages in Fig. |2l The dots mark the ranges of the 
mixed phases. For the PNS matter at low density before the 
onset of quarks, the EOS becomes more and more stiffer with 
the decrease of lepton fraction and the increase of entropy den- 
sity. This mainly attributes to the decreased Yp/Y„, as shown in 
Fig.[T] which excites more neutrons to higher energy states. In 
contrast, the EOS of the mixed phase with a larger lepton frac- 
tion is much stiffer This is because the pure quark phase has a 
stiffer EOS, but the corresponding hadron phase with the same 
lepton fraction and entropy density has a softer one. To fulfill 
the Gibbs condition of the chemical and mechanical equilib- 
rium, the phase transition can only take place at relatively larger 
energy density. For the case with a stiffer hadronic EOS at lower 



3 



density, the Gibbs condition can be realized at lower energy 
density to drive the hadron-quark phase transition. The similar 
results have been obtained when the NJL model is taken 134 1. 
The cold neutron star matter at lower density has an EOS be- 
tween the initial conditions and the end of heating stage, and a 
softest one for the mixed phase after quarks appearing. 



Comparing the results obtained in the NJL model in Il34ll . 
we find the pressure of the mixed phase of a PNS given by 
the PNJL model is larger than that of NJL model. This re- 
flects that the confinement effect (gluon field) is important at 
finite temperature. With the PNJL model, the pressure of quark 
matter at finite temperature are much smaller than that of NJL 
model, therefore the phase transition can only take place at rel- 
atively larger density to assure the quark-phase pressure can 
match that of the hadron phase. The details can be found in 
refs. 1I28, 29] where the hadron-quark phase transition for sym- 
metric and asymmetric matter related to heavy-ion collision ex- 
periment has been investigated in the Hadron-(P)NJL model, 
and the same results were obtained. 

The above illustration also explains the inverse of pressure 
of the mixed phase in the cases of (5 = 2, Yy^ = 0) and 
(S = 0, Yy^ = 0) as given by the NJL model in Q. All the fea- 
tures of equations of state will be reflected in the mass-radius 
relations. 
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Figure 3: Mass-radius relations of PNS and NS at several snapshots along a 
PNS evolution. The solid curves are the results without quarks and the dash 
curves are the results with hadron-quark phase transition. 



Finally, in Fig.|3]we plot the mass-radius relations of several 
snapshots taken above along the evolution a PNS. This figure 
gives us a more intuitive picture about the PNS evolution with 
a hadron-quark phase transition. Firstly, in the heating stage 
the PNS expands with the decrease of lepton fraction (Flp) and 
the increase of entropy density (S). Simultaneously, the center 
density decreases with the expansion of the star. When neutri- 
nos are free, the star begins cooling, and then the star shrinks, 
leading to the center density increasing again. In this stage for 
a star with quarks, a part of nuclear matter transforms to quark 
matter and the EOS becomes softer in the mixed phase. The 
star may collapse into a black hole in the cooling process if the 



EOS in the core cannot resist the gravity. This point is different 
from the result derived in the NJL model. 

The radio timing observations of the binary millisecond pul- 
sar J1614-2230 with a strong general relativistic Shapiro delay 
signature, implies that the pulsar mass is L97 + 0.04 M© I40ll . 
The discovery of this massive pulsar rules out many soft equa- 
tions of state. With the improved quark model, our calcula- 
tion shows the maximum mass of cold neutron star with decon- 
fined quarks is slightly larger than 2 Mq. Because the quark 
model parameters are determined by experiments and lattice 
QCD simulation, a stifFer hadronic EOS is needed to fulfill this 
constraint. On the other hand. If a vector interaction is included 
for quark matter, the EOS will be stiffer and the maximum mass 
of hybrid star will be improved. 

In summary, we have studied the evolution of PNS with a 
hadron-quark phase transition with a more reliable quark model 
including both chiral dynamics and (de)confinement effect. The 
calculation shows the quark phase may exist in the whole pro- 
cess of a PNS evolution. The trapped neutrinos affect greatly 
the ration of Yp / Y„ and the EOS at low density before quarks 
appearing. In the heating stage, with the deleptonization and 
the increase of entropy density, the PNS expands and the cor- 
responding center density decreases. In contrast, in the cooling 
stage, the PNS shrinks and the center density increases. During 
this process a part of nuclear matter transforms to quark matter, 
and the PNS may collapse into a black hole if the EOS in the 
core is not stiff enough. 

In this study, only the qq interaction is considered for quark 
phase because the relevant model parameters can be fixed with 
experiments and lattice QCD simulation. In analogy to BCS 
theory, color superconductivity in low-temperature and high- 
density QCD matter may appear and there may exist rich phase 
diagram. The coupling constant of this interaction channel af- 
fects the equation of state of quark matter and the onset den- 
sity of quark phase in compact star One question is that the 
coupling parameter cannot be fixed from heavy-ion experiment 
or lattice QCD simulation. On the other hand, the (isoscalar) 
vector interaction channel has been included in some studies. 
Such an interaction reduces the effective quark chemical po- 
tential but contributes to the pressure of quark matter This 
interaction stiffens the EOS of quark matter and increases the 
maximum mass of a hybrid star. Compared with the hadron 
Walecka model, the (isoscalar) vector interaction in quark mat- 
ter plays the role corresponding to the u meson. One drawback 
is that the coupling constant is usually taken as a free parameter 
and its strength affects greatly the Critical End Point of chi- 
ral symmetry restoration relevant to heavy-ion collision exper- 
iment. Besides, the isovector vector interaction (corresponding 
to the role of p meson in hadron phase and influencing the sym- 
metry energy of quark matter) can be also introduced for asym- 
metric quark matter. For these interaction channels as well as 
hyperon degrees of freedom, the difficulty is the uncertainties 
of the relevant coupling parameters with the lack of experiment 
data, so we temporarily omit these interactions in this study. 
And a systematic investigation on these problems is in progress 
as a further study. 
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